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NOTATION 


leading-edge thickness 

leading-edge thickness defined by equation (11) (see sketch (a)) 
leading-edge thickness defined by equation (12) (see sketch (a)) 
variable contained in equation (Al) 
free-stream Mach number 

bluntness reduced Mach number (determined from the ratio of total 
pressure behind the shock at the leading edge and the free- 
stream static pressure) 

bluntness Reynolds number, (R/inchJ^b) 

transition Reynolds number, (R/inch) oo (x-t) 

transition Reynolds number for b = 0 as defined by equation (1) 
transition Reynolds number as defined by equation (3) ■ 


transition Reynolds number which accounts for both bluntness 
reduced Mach number and reduced unit Reynolds number effects, 
(R/ inch) oo 




(R/inch) n 


transition Reynolds number which accounts for only bluntness 
reduced unit Reynolds number effects, (Rt,o) (R^nch)*” 


free-stream Reynolds number per inch, — ^ 


free-stream static temperature 
total temperature 

static temperature based on bluntness reduced Mach number (M n ) , 
To 

1 + 0,2 M n 2 
wall temperature 
free-stream velocity 

velocity based on bluntness reduced Mach number (M n ) , (49 /T n )M n 
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distance from leading edge to end of transition 

Y c height above surface where maximum fluctuation energy is indicated as 
defined in reference 1 (see chart 4 herein) 
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ratio of critical layer to boundary- layer thickness used in defining 


b l > 


essentially. 


Y c 6? 

-r — = — (see chart 4 and eqs. 

°t 6 t 


(A7) and (A8)) 


thickness of bluntness reduced Mach number and unit Reynolds number 
layer (see sketches (b) and (c)) 


Ysb perpendicular distance to the horizontal plane of symmetry from the 
sonic point on the body (see sketches (b) and (c)) 


y ratio of specific heats (assumed to be 1.4 for air) 

boundary-layer thickness 
boundary- layer displacement thickness 
A leading-edge sweep 

u viscosity based on T 

y n viscosity based on T n 

v n kinematic viscosity based on T n 
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CHARTS FOR ESTIMATING BOUNDARY- LAYER TRANSITION 


ON FLAT PLATES 

Edward J. Hopkins, Don W. Jillie, and Virginia L. Sorensen 

Ames Research Center 


SUMMARY 


Charts are presented for rapidly estimating the end of boundary- layer 
transition for flat-plate wind-tunnel models with supersonic leading edges at 
an angle of attack of 0°. The charts were developed from the semiempirical 
method of Deem and Murphy who derived an equation that accounts for the com- 
bined effects of Mach number, unit Reynolds number, leading-edge sweep, 
leading-edge bluntness, and wall temperature on transition Reynolds number. 


INTRODUCTION 


For wind-tunnel tests, it is necessary to know the condition of the 
boundary layer so that the results can be interpreted correctly for extrapola- 
tion to flight Reynolds numbers. In addition, at hypersonic Mach numbers it 
is important to know whether the boundary layer is laminar or turbulent both 
from the standpoint of drag and of heat transfer. 

Until recently no integrated method existed for estimating the end of 
boundary- layer transition since so many factors were known to affect transi- 
tion; therefore, emphasis in past investigations was placed on a separate 
evaluation of each factor. Many of these factors are discussed in references 
2 and 3. Deem and Murphy developed a semiempirical equation in reference 2 by 
which five important variables known to affect transition can be taken into 
account: Mach number, unit Reynolds number, leading-edge sweep, leading-edge 

bluntness, and wall temperature. The method was developed from data on 
"aerodynamical ly" smooth models having surface roughness less than that 
required to influence boundary- layer transition. Since the data were taken in 
wind tunnels with different turbulence levels, the results should be consid- 
ered as only representative average transition Reynolds numbers. The method 
is, furthermore, only applicable to flat -plate models with supersonic leading 
edges at an angle of attack of 0°. Examination of 291 experimental points in 
reference 2 gave a standard deviation of 33 percent between the measured and 
calculated transition Reynolds numbers. However, much of this deviation can 
be explained on the basis of variations in wind-tunnel turbulence and impre- 
cise measurements of transition and leading-edge bluntness. 

The charts presented are based on the semiempirical method of reference 2 
and offer a simple method for determining the approximate end of transition 
for flat plates mounted in wind tunnels. The charts cover ranges of 
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leading-edge sweep from 0° to 80°, unit Reynolds numbers from 10 3 to 10 7 per 
inch, Mach numbers from 1.1 to 12, and leading-edge thicknesses from 0 to 
0.5 inch. 

Although the charts presented are based on wind-tunnel transition results, 
they should be useful in making a first approximation of the minimum amount of 
laminar flow to be expected in flight. In general, transition Reynolds numbers 
in flight are expected to be higher than in wind tunnels, as evidenced by fig- 
ure 3 of reference 4. By firing a hollow cylinder into still air, James (ref. 
5) obtained nearly four times the transition Reynolds number as did the present 
authors for a flat plate mounted in a wind tunnel, although both models had 
nearly the same leading-edge thickness. Moreover, transition Reynolds numbers 
as high as 33 million were measured at a Mach number of 3.15 when large highly 
stabilized cone-cylinder rockets were fired into the atmosphere (ref. 6) . 


DESCRIPTION OF THE DEEM AND MURPHY METHOD 


The important concepts and assumptions involved in the Deem and Murphy 
method for predicting transition Reynolds number will be described briefly. 

Some sketches will be included to help in the interpretation of certain blunt- 
ness criteria. The method is restricted to flat-plate wind-tunnel models hav- 
ing aerodynamic ally smooth surfaces. In addition, the method as developed is 
only applicable to models at an angle of attack of 0° with supersonic leading 
edges of nearly semicircular cross section. The method should be applicable, 
however, to models at small angles of attack provided the local flow conditions 
above the model are used in the equations. (If the leading edge is subsonic, 
the transition process is more complex because of the formation of a leading- 
edge "bubble," see ref. 4.) The method accounts for the effects of five quan- 
tities known to effect transition: Mach number, leading-edge bluntness, 

leading-edge sweep, wall temperature, and unit Reynolds number. 


Mach Number 

Whitfield and Potter (ref. 7) showed that transition Reynolds number for 
zero leading-edge bluntness (Rt,o) increased over five times when the Mach 
number was increased from about 3 to 8. This increase in Rt,o cannot be 
explained on the basis of the physical leading-edge bluntness reducing the unit 
Reynolds number (Moeckel effect, ref. 8), since the values of Rt,o were 
obtained by extrapolation of data curves to zero bluntness. Deem and Murphy 
(ref. 2) empirically derived an equation from existing data, after reducing the 
data to a common unit Reynolds number of 3xl0 5 per inch, to give this Mach 
number effect as 


R t ,o = l x 10 6 + 0.36xl0 6 |M ot - 3 | 3/2 (1) 

Whitfield and Potter found that Rt did not reach a limit even though 
some of their data for models with blunt leading edges were obtained close to 
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the bluntness Reynolds number (Rb) required according to Moeckel f s analysis 
for maximum influence (i.e., Rb = 3000). To account for this insufficiency of 
the Moeckel effect alone giving the correct prediction of Rt, Whitfield and 
Potter empirically derived the following approximate formula for transition 
Reynolds number 


R t ,i = ( R t,o for MJ 


(R/inch)^ /R t ,o for M n ' 
(R/inch) n \ R t,o for M a 


( 2 ) 


Deem and Murphy’s 
number becomes 


( R t,o 


equation (1) written in terms of the bluntness reduced Mach 
for M n ) = (Rt,o) n = lxlO 6 + 0.36xl0 6 |M n - 3| 3/2 (3) 


In the Deem and Murphy method, it is assumed that equation (2) applies 
for relatively sharp leading edges, but that for very blunt leading edges only 
the reduced unit Reynolds number (Moeckel effect) applies; therefore, the 
equation for R t becomes 


(R/ inch)^ 

R t,2 = ( R t,o for MJ - (4) 

(R/inch) n 


in which (R^ 0 for MJ = Rt,o from equation (1). In the next section the 
leading-edge bluntness criteria for using equation (2) or (4) will be 
discussed. 


Leading-Edge Bluntness 

A typical curve representing the variation of Rt with bluntness for an 
unswept plate having a semicircular leading edge was derived by Deem and 
Murphy from a limited amount of data and is shown in sketch (a) . Deem and 



Sketch (a) 


Murphy assumed that Rt,l was pre- 
dictable by the Wl.itfield-Potter 
equation (2) and that Rt,2 was pre- 
dictable by the Moeckel equation (4) . 
The value of Rt,o is determined by 
the empirical equation (1) . To find 
the transition Reynolds number within 
regions 1 and 2 of sketch (a) it was 
assumed that a linear interpolation 
could be made as shown. It follows 
that the equation defining the tran- 
sition Reynolds number for region 1 
is 
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for region 2 


R t = R t ,0 + ( R t,l - R t,o) 


C5) 


R t = R t,i 



C R t,2 - R t,l) 


( 6 ) 


and for region 3, 


(R/inch)^ 

R t = R t,o 

(R/inch) n 


(7) 


In order to use equations (5) through (7) it is necessary to establish 
the bluntness criteria Ct> i and ^> 2 ) required to define the regions of sketch 
(a) . Deem and Murphy accomplished this by analyzing a limited amount of exper- 
imental data. They found that good correlation is obtained for predicting 
transition Reynolds number provided b]_ and b 2 are defined in terms of the 
thickness of the bluntness reduced Mach number layer (Y n ) relative to the 
boundary- layer thickness. For bj , the displacement thickness was accounted 
for in defining the critical thickness of the reduced Mach number layer, but 
for b 2 the displacement thickness was ignored. Geometrical details 1 and 
equations defining b^ and l >2 are presented with sketches (b) and (c) , 
respectively. 



Sketch (b).- Geometrical details for bi« 


1 In either sketch (b) or (c) , M n is the bluntness reduced Mach number 
determined from the ratio of the total pressure behind the normal shock at the 
leading edge to the free-stream static pressure. After Moeckel, M n is con- 
sidered to be nearly the same as the Mach numbers within the shear- layer height 
CYn) , and, therefore, is used in this manner to simplify the calculations. 
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Sketch (c).- Geometrical details for bp. 


From sketch (b) , following Moeckel's assumption regarding the critical 
bluntness for full bluntness effect on transition, for a flat-plate model with 
a semicircular leading edge 2 * * * 

Y SB = | sin 51.8° (8) 

Y n = " 6 t O) 


Dividing equation (9) by (8) and solving for b, we obtain 


b 




( 10 ) 


Using experimental data and equation (10) for which bluntness was a variable 
and Y n /YsB was considered a parameter. Deem and Murphy found that bi could 
be determined provided Y n /Y 5 B = 3, a value found to be independent of free- 
stream Mach number. Therefore, the equation for b^ becomes 


bi 




(ii) 


2 The development in reference 2 is based on an assumed semicircular 

leading edge. Such a leading edge is generally representative of a so-called 

sharp leading edge and small departures from such a shape will not alter the 

results greatly. 
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This bluntness criterion differs from that of Moeckel (ref. 8) for which 
Y n /YsB was shown to vary with free-stream Mach number. This difference is 
related, evidently, to Deem and Murphy's choice of iterating on equation (10) 
with experimental data and by obtaining best correlation by using two differ- 
ent bluntness criteria, t>! and b 2 - 

For models with very blunt leading edges relative to the transition-point 
boundary- layer thickness (region 3 of sketch (a)), the experimental transition 
data examined by Deem and Murphy indicated good correlation if the displace- 
ment thickness is ignored in the definition of b 2 . For this case, as shown 
in sketch (c) , the reduced Mach number inviscid shear layer is measured from 
the model surface. The ratio Y n /YsB was again considered to be a parameter 
but was found, from experimental data and equation (10) , to have a value of 
1.0, independent of Mach number. According to Deem and Murphy's empirical 
findings the equation for b 2 can be written from equation (10) (with 6^ 0 

and Y n /YsB = 1) as 


bo = 


L 2 

(sin 51 . 8° / \ Y n /Y SB ) 


= 2.545 6* 


( 12 ) 


Leading-Edge Sweep 

Leading-edge sweep can affect boundary- layer transition in two ways: 

First, when the leading edges are supersonic, the Moeckel effect (reduced unit 
Reynolds number due to leading-edge bluntness) becomes progressively smaller 
as sweep is increased because of the reduced strength of the leading-edge 
shock; that is, a constant assumed transition Reynolds number based on local 
flow contributions corresponds to a progressively smaller transition Reynolds 
number based on free-stream conditions as sweep increases, because the local 
unit Reynolds number also increases. The validity of the Moeckel assumption 
as applied to sweep was demonstrated in reference 4 in which it was shown that 
when the leading edge is supersonic, the transition Reynolds number decreases, 
at least qualitatively, with increases in sweep, as predicted, and that at a 
Mach number of 0.27 the transition Reynolds number remains nearly constant 
with increases in sweep to about 30°. Second, with both subsonic and super- 
sonic leading edges, the variable crossflow within the boundary layer creates 
a twisted boundary- layer profile that can lead to boundary- layer instability 
and transition when a certain critical crossflow Reynolds number is reached, 
as suggested by Owen and Randall (ref. 9) . Evidence that this phenomenon does 
occur is given in references 9 and 10. 

Deem and Murphy semiempirically accounted for both the above sweep effects 
in the following manner: 

1. In all t hree regions of sketch (a), Rt,o or (Rt,o) n is multiplied by 
the factor /cos A. The reason this factor improved the correlation of the data 
is probably related to a decrease in effective bluntness as sweep increases. 
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2. In regions 2 and 3 at sweep angles greater than 25°, it is assumed 
that crossflow is dominant in affecting transition and that the total ef fect 
of sweep on transition can be approximated by multiplying Rt,o by the /cos A 
factor without any further corrections for bluntness on either Mach number or 
unit Reynolds number. 

3. In region 2 for sweep angles less than 25° and in region 1 for all 
sweep angl es, both crossflow and bluntness are considered important; therefore, 
the /cos A factor is used with the bluntness reduced Mach number and unit 
Reynolds number factors. 

4. In region 3 for sweep angles less than 25°, only the /cos A factor is 
used with the bluntness reduced unit Reynolds number. 


Wall Temperature 

The effect of wall temperature is accounted for in the method indirectly 
through the boundary- layer thickness equation used to define the bluntness 
criteria. Thus, by this treatment the transition Reynolds number is affected 
by wall temperature only in regions 1 and 2 and not at all in region 3 (with 
b >> b 2 ) . 


Unit Reynolds Number 

According to an analysis of existing data made by James (ref. 5) , the 
variation of transition Reynolds number with unit Reynolds number can be 
expressed as 


loglO R t = Ci + 0.4 logic (R/inch)^ (13) 

where Ci is dependent on many variables. In the Deem and Murphy method Ci 
is assumed to depend on the four variables: Mach number, leading-edge blunt- 

ness, leading-edge sweep, and temperature; and the second term in equation (13) 
is used to predict the effect of unit Reynolds number. 

For completeness, the equations used for the calculations are summarized 
in appendix A. A more detailed development of the equations is given in 
reference 2. 


PRESENTATION OF RESULTS 


Chart 1 gives the effect of Mach number (M = 1.1 to 12) on transition 
Reynolds number for various leading-edge thicknesses (b = 0.0001 to 0.5 inch) 
and unit Reynolds numbers ((R/inch) = 10 3 to 10 7 ) . Chart 2 gives the effect of 
Mach number (M^ = 1.1 12) on the normalized transition Reynolds number for 

various angles of sweep (A = 0 to 80°), unit Reynolds numbers 
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((R/inch)^ = 10 3 to 10 7 ) and leading-edge thicknesses (b = 0.0001 to 0.5 inch.) 3 * 4 
In both sets of charts the wall temperature and the total temperature were held 
constant at 400° R and 500° R, respectively. A limited study of the effects of 
large temperature changes on the transition measured on flat plates indicated 
that temperature effects were small and were adequately accounted for by the 
slight indirect effect on boundary- layer thickness as contained in the method. 
This may seem contrary to the commonly accepted large variations in transition 
due to surface cooling; however, the latter effects may be related to the three 
dimensionality of the flow. Calculated transition results presented in chart 3 
for two different Mach numbers and unit Reynolds numbers for the unswept case 
show a very small effect due to temperature changes. Since the method should 
be considered only for obtaining approximate values of transition Reynolds num- 
ber, no additional corrections for this secondary effect of wall-temperature 
ratio appear justified. 5 


Sample Calculation 


It will be assumed that the following conditions exist on a flat plate: 

M =6.0 

oo 

(R/inch) ro = 10 6 

b = 0.01 inch 
A = 60° 

From chart 1(d), (Rt)^_ 0 = 1.4xl0 7 and from chart 2(d), Rt/( R t)A- 0 = 0.23 
for the above conditions. The transition Reynolds number is 


3 Note that for a particular sweep angle the minimum Mach number shown is 
that at which the leading edge first becomes sonic. It is shown in reference 4 
that when the leading edge becomes subsonic, transition moves close to the 
leading edge because of the influence of a localized leading-edge separation 
"bubble." 

4 The sharp discontinuities in some of the curves presented in charts 2(d) 
and 2(e) are related to the arbitrariness in defining the bluntness, b. In 
real flow, the curves would probably have a smoother variation with Mach number 
than shown. 

5 Chart 2 should be applicable to flight vehicles as well as wind-tunnel 
models. Chart 1, however, should be restricted to wind-tunnel models or used 

to estimate the minimum expected transition Reynolds number for flight, because 

the transition Reynolds numbers given were derived from wind-tunnel experiments 
in which wall disturbances and free-stream turbulence were present. The 
usually large difference in static temperature for wind-tunnel and flight 
models, however, is not expected to have an important influence on transition, 
provided the surfaces are relatively flat. 
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R t - Ma=o [ R t/ C^t) a =0 ] - 3.2*106 


CONCLUDING REMARKS 


From the charts presented, a rapid estimate can be made of the boundary- 
layer transition on flat-plate models mounted in wind tunnels. These charts 
are restricted to models with supersonic leading edges at an angle of attack 
of 0°. Mach number, unit Reynolds number, leading-edge bluntness, sweep, and 
wall temperature are the variables considered. The wall temperature was shown 
for flat plates to be of secondary importance in affecting transition; there- 
fore, the charts are presented for a single wall temperature. 

It should be emphasized that these charts give only an approximate 
transition Reynolds number for flat plates mounted in wind tunnels. Since in 
flight the maximum obtainable transition Reynolds numbers should be consider- 
ably higher than those for a wind tunnel, the charts should be used only to 
obtain a first approximation of the minimum transition Reynolds number 
expected for flight vehicles. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, April 2, 1970 
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APPENDIX A 


EQUATIONS USED FOR CALCULATING TRANSITION REYNOLDS 
NUMBER BY THE DEEM AND MURPHY METHOD 

BASIC TRANSITION REYNOLDS NUMBER EQUATION 


The basic equation for calculating transition Reynolds number is 


logic R t = c i + 0.4 log 10 (R/inch)^ 


(Al) 


where Ci depends on Mach number, leading-edge bluntness, leading-edge sweep, 
and wall temperature. The second term in equation (Al) is considered indepen- 
dent of Ci and accounts for the unit Reynolds number effect. 


REDUCED MACH NUMBER AND UNIT REYNOLDS NUMBER EQUATIONS 


In the method, transition Reynolds number is affected by changes in Mach 
number and unit Reynolds number near the surface resulting from leading-edge 
bluntness and the associated leading-edge shock losses. In reference 4, 
expressions were given for these reduced numbers for which y = 1.4 was 
assumed. 1 The same expressions are given below: 


1/2 


M n = 


5/2 


l 7M 2 cos 2 A 


6M 


cos 2 A(M 2 + 5) 


5 (M 2 cos 2 A + 5) 


7/2 


2/7 


(A2) 


and 


where 

y 

y 


(R/inch) j 

(R/inch) c 

0.03665 T (T n + 198.6) 


1/2 


M 

_r 

M 


J n 


n 


(T n ) 


3/2 


for T < 200° R and T n > 200° R 


(A3) 


(A4) 


x No explicit expressions for reduced Mach number or unit Reynolds number 
can be written when real gas effects are taken into account. Since the method 
of reference 2 is highly approximate and the real gas effects are known to be 
small, ignoring these effects seems justified. 
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It is necessary to establish the value of the leading-edge bluntness, b, 
relative to the bluntness criteria, bj and b 2 , defined by equations (11) and 
(12). If Creager's equation is used for boundary- layer thickness as given in 
reference 11 and Y n /Ygg = 3.0 (as determined in ref. 2), equation (11) for 
small bluntness becomes 



for T w > 200° R and <; 200° R 


(A8) 2 » 3 

2 Since the unknown transition Reynolds number, R-^ occurs in equations 
(A7) , (A8) , (A10) , and (All), it is necessary to obtain R t by iteration with 
equation (Al) which contains . In turn, the proper equation for C^, given 
hereinafter, can only be determined after bj and b 2 have been calculated. 

3 In equations (A8) and (All) the linear approximation to Keyes' viscosity 
equation (ref. 12) was used for temperatures equal to or below 200° R; whereas 
Sutherland's viscosity formula (ref. 13) was used for temperatures above 
200 ° R. 
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where Y c /<5 t , as noted by Potter and Whitfield (ref. 1), is, essentially. 



Values of Y c /6-£ from reference 1 are presented in chart 4. If Creager's 
equation is modified for boundary- layer thickness to include the effect of 
bluntness on the local surface conditions (M m M n , T n , and 

(R/inch)^ (R/inch) n ) and if Y n /YsB = 1.0 (as determined in ref. 2) is used, 
equation 00 (12) for large bluntness 4 becomes 



or 

b 2 = 2.545 ^ + 0.1328 + (M^ 2 ) 

\Mn 2 T n M n 2 ) 


for T w > 200° R and T n < 200° R 

(" (R/inch)! 

(R/inch)^ |j R - /inch )J (All) 2 » 3 

TERM C 2 FOR EQUATION (Al) 

Since the transition Reynolds number given by equation (Al) is dependent 
on Ci which depends on sweep and the bluntness criterion (b relative to 
b 2 and b 2 ) which also contains Rp, it is necessary to use an iterative pro- 
cess to calculate Rp. The following equation defines Ci and the required 
sweep and bluntness criterion: 

2See footnote, p. 11. 

3 See footnote, p. 11. 

4 For this case, the entire boundary layer is engulfed in the Mach number 
reduced inviscid shear layer in which the Mach number is nearly constant and 
equal to M n . 
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For 0 < b - bi, all sweep angles 


Cl = log 10 |(lxl0 6 + 0.36x10 s Im^ - 3| 3/2 )(cos A) 1/2 


(R/inch) 


1 + 


00 / b 


QxlO 6 + 0 . 36xl0 6 | M n - 3| 3/2> 


(R/mch) n \biA lxl0 6 + o.36xl0 6 IM^ - 3| 3/2 / bl 


i- • 


19 


For bl < b < b 2 , o < A < 25 c 


C 1 - logio 


(lxlO 6 + 0.36x10 s |M n - 3| 3 / 2 )(cos A) 1/2 


(R/inch) OT ^ ( 

inch) n 1 + \t 


b - bA/lxlO 6 + 0 . 36x10 s |m - 3 | 3/2 \ / b - b x N 


(R/inch) 


2 bl /Vlxl0 s + 0.36x10 s |M n - 3 | 3 / 2 / \b 2 - 


For b > b 


C 1 = l°gl0 


2 j 0 < A < 25° 

I (l x 1 0 6 + 0.36x10 s |M m - 3 | 3/ 2 ) (cos 


(R/inch) ^ 
(R/inch) n 


2.19 


For b > bi, A > 25° 

Ci = logi 0 |^(lxl0 6 + 0.36x10 s IM^ - 3 1 3/ 2 ) (cos A) 1/2 j - 2.19 


(A12) 


- 2.19 
(A13) 

(A14) 

(A15) 
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Chart 1. - Estimated effect of Mach number on transition Reynolds number 
for a flat plate with no leading-edge sweep; T w = 400° R, T 0 = 500° R. 


15 











Lart Es 
number fo: 
T w = 400° 


4 



5 


10 

0 . 00 < 

o.oo; 

o.oi 



















Chart 2.- Concluded. 
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(a) (R/inch)^ = 10 3 


Chart 3 .- Estimated effect of wall -to -total temperature ratio on transition 
Reynolds number for a flat plate; b = 0.01 inch, A = 0, T w = 400° R 
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